Disentangling glass and jamming physics in the rheology of soft materials 
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The shear rheology of soft particles systems becomes complex at large density because crowding ef- 
fects may induce a glass transition for Brownian particles, or a jamming transition for non-Brownian 
systems. Here we successfully explore the hypothesis that the shear stress contributions from glass 
and jamming physics are 'additive'. We show that the experimental flow curves measured in a large 
variety of soft materials (colloidal hard spheres, microgel suspensions, emulsions, aqueous foams) 
as well as numerical flow curves obtained for soft repulsive particles in both thermal and athermal 
limits are well described by a simple model assuming that glass and jamming rheologies contribute 
linearly to the shear stress, provided that the relevant scales for time and stress are correctly iden- 
tified in both sectors. Our analysis confirms that the dynamics of colloidal hard spheres is uniquely 
controlled by glass physics while aqueous foams are only sensitive to jamming effects. We show 
that for micron-sized emulsions both contributions are needed to successfully account for the flow 
curves, which reveal distinct signatures of both phenomena. Finally, for two systems of soft mi- 
crogel particles we show that the flow curves are representative of the glass transition of colloidal 
systems, and deduce that microgel particles are not well suited to studying the jamming transition 
experimentally. 



I. INTRODUCTION 

The emergence of solidity in disordered assemblies of 
soft particles is observed in large variety of systems, ev- 
eryday examples including toothpaste, shaving foam, or 
paints [H, 0- Because of their practical use, a fundamen- 
tal understanding of the flow property of dense suspen- 
sions is necessary, but this continues to represent a con- 
siderable challenge to physicists. In particular, achieving 
a detailed understanding of fluid-to-solid transitions in 
amorphous materials is recognized as an important chal- 
lenge for modern condensed matter physics [3, Q • 

Two types of dense particulate systems have been con- 
sidered. One class is composed of particles that are small 
enough to be sensitive to thermal fluctuations. Brownian 
forces are for instance relevant in colloidal suspensions of 
micron-sized particles or smaller. In that case, colloids 
undergo collisions with the solvent molecules which play 
the role of a thermal bath. The emergence of amorphous 
solids at large density in such materials is usually called 
colloidal glass transition Q , by analogy with the molecu- 
lar glass transition observed in organic or polymeric liq- 
uids 0. 

In a second class of systems thermal fluctuations play 
a negligible role, with aqueous foams being one exam- 
ple. In these non-Brownian systems, the particle size is 
typically very large (say, in the millimeter range), such 
that the effect of collisions with solvent molecules is so 
weak that it can be safely neglected on typical experi- 
mental timescales. The emergence of solidity in athermal 
amorphous systems was termed jamming transition Q, 
by analogy with dry granular media Q ■ 

Surprisingly, despite large differences in the underlying 
microscopic dynamics, the phenomenology of glass and 
jamming transitions appears remarkably similar, inpar- 
ticular when soft particle systems are considered Q . 



Below some critical density, the system can flow un- 
der infinitesimal shear stress. The nature of the flow 
is Newtonian, and the system thus behaves as a viscous 
fluid. Above that critical density, the system does not 
flow if the shear stress is below a characteristic finite 
value, called the yield stress. The system is now a soft 
solid. Close to the critical density, the rheology is typ- 
ically strongly non-linear: it is very sensitive to small 
density changes, and the microscopic dynamics becomes 
spatially and temporally heterogeneous [l(| ■ 

This apparent similarity suggests that theoretical 
progress for one type of systems might impact research 
for the other. For Brownian systems, flow curves are 
typically analyzed by assuming that the glassy dynam- 
ics arising in thermal suspensions at rest is disrupted by 
the externally imposed shear flow flll - fTij . A number of 
experimental studies, in particular in the colloidal liter- 
ature, have then used these ideas to organize and quan- 
titatively describe flow curves measured in the vicinity 
of the glass transition 15|. By contrast, much less is un- 



derstood regarding the microscopic dynamics of athermal 
suspensions near the jamming transition, because the im- 
posed shear flow represents at the same time the external 
driving force and the source of microscopic fluctuations 
responsible for particle motion. Nevertheless, detailed 
numerical analysis of model systems and theoretical mod- 
elling have suggested specific functional forms and scaling 
procedures to or gan ize rheological data around the jam- 
ming transition [Ia - [l9j . In fact these ideas have been 
used in a larger variety of systems, including atomistic 
glasses [2(| and colloidal particles [2l|, where their rele- 
vance is not obvious a priori. 

Recently, we have used numerical simulations to 
study the rheology of dense assemblies of soft repulsive 
spheres [22[. Despite its simplicity, this model is use- 
ful as it is known to display glassy dynamics at finite 
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temperatures [23| . and to undergo a geometric jamming 
transition at zero temperature [8|, [2_^. By varying the 
relative strength of energy dissipation and thermal ag- 
itation over a broad range, we have demonstrated that 
glass and jamming physics impact the steady state flow 
curves over distinct stress scales and time scales, which 
can be varied independently. These results confirm that 
the two phenomena are actually distinct, and that their 
respective influences on rhcological behaviour should not 
be confused. In addition, we have suggested that the 
numerical flow curves are well described by a simple 'ad- 
ditive' model, where the glass and jamming contributions 
to the shear stress can be captured separately and then 
combined linearly into a unified model [22I ] . 

In the present work, we expand our presentation of 
this additive model to explain in greater detail its con- 
struction and application to our numerical data. We then 
confront our rhcological model deduced from the analy- 
sis of a simplistic soft particle system to a broad range 
of experimental data taken from the literature. Overall, 
our analysis validates the hypothesis that, despite their 
strongly nonlinear nature, glass and jamming physics lin- 
early contribute to the rhcology of soft particle systems 
and act over distinct sectors, whose relative importance 
depends on the experimental system at hand. We also 
show that our model and theoretical analysis provide use- 
ful tools to organize and explain rheological data stem- 
ming from different experimental sources and materials. 
These tools can be used to efficiently identify which of 
the glass or jamming contributions is most relevant, or 
whether both types of physics are in fact needed to cor- 
rectly describe the data. 

Our paper is organized as follows. In Sec. [TT] we con- 
struct our additive rheological model which we apply to 
numerical data obtained from simulations of harmonic 
spheres. In Sec. IIIII we use this model to analyze steady 
state flow curves measured in colloidal hard spheres, mi- 
crogel suspensions, emulsions, and aqueous foams. In 
Sec. lIVi we summarize our findings and discuss some per- 
spectives for future work. 



II. 'ADDITIVE' MODEL FOR GLASS AND 
JAMMING RHEOLOGY 

A. Flow curves for harmonic spheres 

To construct a general model for the complex rheol- 
ogy of dense suspensions, we use the numerical results 
obtained for a specific model as a reference. In previ- 
ous work, we have studied the steady state rheology of 
harmonic spheres in the overdamped limit [22[ . This can 
be seen as a simple model to describe the physical be- 
havior of dense suspensions of deformable spheres, such 
as the emulsions and colloidal suspensions considered in 

sec. urn 

Since the details of the model were already reported 
in our previous work, we only briefly summarize the key 



ingredients. We consider N harmonic spheres contained 
in a volume V. We use Lees-Edwards periodic bound- 
ary conditions [25| and solve the following equations of 
motion numerically: 



dt 



ivie x ) 



dv{\fj - fj\) 

3=1 



(1) 



Here rj represents the position of particle i, its y- 
componcnt, and e x the unit vector along the x-axis. The 
damping coefficient, £, and the random force, Ri(t), obey 
the fluctuation dissipation relation: (Ri ja (s)Rj t p(s')) = 
2kBT£5ij5 a pS(s — s'), where T is the temperature and 
ks is Boltzmann's constant. The interaction potential is 
a purely repulsive harmonic interaction truncated at the 
particle diameter, 



(2) 



where a is the particle diameter and O(x) is the Heaviside 
function. To avoid crystallization issues, we work with a 
50:50 binary mixture of particles with diameter ratio 1.4, 
but this is otherwise largely irrelevant. 

By construction, the system has two characteristic en- 
ergy scales, namely the thermal energy, fcgT, and the in- 
teraction energy of particles e. The ratio of these energy 
scales, kgT/e, is an important control parameter. Phys- 
ically, it is a measure of the particle softness, expressed 
in units of kinetic energy. 

Glass and jamming effects can be distinguished most 
clearly in the low-softness limit fcgT/e —> 0. This can be 
obtained in two ways. One can either take the limit of 
vanishing kinetic energy, T — ¥ 0, at fixed repulsion en- 
ergy scale e. One then studies the jamming transition of 
athermal packings of soft repulsive spheres. The alterna- 
tive is to keep temperature constant and send e — > 00, a 
limiting case which corresponds to studying the physics 
of the thermalizcd hard sphere fluid. 

The two energy scales fc^T and e naturally provide two 
characteristic time scales and stress scales. The micro- 
scopic time scale for Brownian motion is 



T 



(3) 



and represents the time it takes a particle to diffuse by 
Brownian motion over a length scale comparable to its 
size. The characteristic time scale for energy dissipation 
is given by 



£a 2 k B T 

t = = TT- 

e e 



(4) 



Likewise, we can define two stress scales: a typical 
stress created by thermal fluctuations, 



<7T 



(5) 
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and an athcrmal stress scale: 



<?0 



-(Tx- 



(6) 



The above expressions make it clear that these time scales 
and stress scales are comparable when particle softness is 
large enough, fc^T sa e, while they become increasingly 
separated as the athermal limit is approached, ksT <C e. 

In terms of the time and stress scales above we can 
now give a more explicit statement of when the limit 
ksT/e — > corresponds to athermal or thermal be- 
haviour. If the limit is taken by lowering T at constant e 
then the shear rate in athcrmal units, 7T0, stays fixed: we 
are exploring the 'athcrmal sector'. Here we then also ex- 
pect stresses to be on the athermal scale, with a/ao of or- 
der unity. If on the other hand we fix T and take e — > 00, 
then the shear rate on the thermal scale, 7T7;, remains 
finite. One then probes the 'thermal sector', with corre- 
spondingly a far = 0(1). Numerically, we find (see be- 
low) that the prefactors in these estimates are such that, 
below and around the glass transition, the interesting 
physics typically takes place for jtt < 1 and <j/<jt < 1- 
As a matter of general orientation we can therefore say 
that phenomena where the shear rate and stress in ther- 
mal units are below unity belong to the thermal sector, 
while the athermal sector corresponds to values of these 
quantities significantly above one. A numerical analysis 
along the same lines has appeared recently (2(| . 

In Fig.[T]we present flow curves obtained by measuring 
the average shear stress, a, under steady state conditions 
created by a constant applied shear rate, 7, and varying 
the packing fraction, <p } and the temperature, T, across 
a broad range. Note that we use stresses and shear rates 
normalized by tt and or as time and stress units; these 
are the appropriate scales for thermal systems. This leads 
us to introduce a 'rcnormalized' shear rate, 



Pr 



7T T , 



(7) 



which is also called the Peclet number. Because the re- 
sults shown in Fig. [I] were discussed in detail in Ref. [22| . 
we only summarize the most salient features of these 
data. 

• High temperature / Soft particles (such as ksT/e = 
10~ 4 ). A single fluid-to-solid transition is ob- 
served in this thermal regime. When the density is 
low, Newtonian flow is observed at low shear rates 
7 <C t^ 1 . With increasing density, the Newto- 
nian viscosity grows and a shear stress plateau ap- 
pears at finite shear rate, corresponding to a strong 
shear-thinning regime. Near the transition density, 
<p « <pg, the viscosity becomes very large, and a 
finite yield stress appears, which grows smoothly 
with density upon further compression. 

• Low temperature / Hard particles (such as ksT/e = 
10~ 7 ). The glass and jamming transition can be 
observed separately When the density is low, 
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FIG. 1: Flow curves of harmonic spheres for temperatures 
from ksT/e = 1CP 4 down to 10 -7 measured in the simula- 
tions first described in Ref. [13] (symbols). In each panel, 
packing fraction increases from bottom to top in the range 
<p = 0.55 — 0.70. The red lines are fits using the additive 
model presented in this work, with parameters listed in Ta- 
ble[TJ The dashed lines delimit the thermal/glass (bottom left) 
and athermal/jamming (top right) sectors. 



two Newtonian regimes appear for 7 <C t t 1 and 
t^ 1 <C 7 <C Tq 1 , respectively. The first Newtonian 
viscosity becomes very large near the glass transi- 
tion density <pc where the system acquires a finite 
yield stress, of the order of the thermal stress scale 
<7t- However the second Newtonian viscosity re- 
mains finite at the glass transition density. On fur- 
ther increasing the density, the second Newtonian 
viscosity then diverges at the jamming transition 
density ipj. At the same time, the yield stress value 
increases rapidly with density; above the jamming 
transition it reaches a value which is controlled by 
the athermal stress scale cto- 

These numerical results directly illustrate that glass 
and jamming transitions represent distinct ways for the 
system to form an amorphous solid phase, since the fluid- 
to-solid transitions occur over different time windows 
(or, equivalcntly, different shear rates) and give rise to 
solids with yield stress having different scales and differ- 
ent physical origins. By tuning the temperature of our 
system, we can move from Brownian soft particles un- 
dergoing a glass transition in the thermal sector of Fig. Q] 
to non-Brownian ones undergoing a geometric jamming 
transition in a different sector. Interestingly, both types 
of physics compete and affect the flow curves for inter- 
mediate values of the temperatures. These observations 
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indicate that the particle softness, ksT/e, is a central 
control parameter for the rheology of soft materials. In 
particular, we note that when the particles are too soft 
(equivalently, when temperature is too large), the jam- 
ming transition has no effect on the rheological data. 



B. Construction of an additive model 

1. The additive hypothesis 

Using the simulation results as a reference, we now con- 
struct a simple rheological model. The major assumption 
is the additivity of the contributions from thermal and 
athcrmal parts to the total shear stress: 



v( T , <P, 7) = 0"G (T, tp, 7) + cr.j (tp, 7) + r) s j. 



(8) 



Here <tg is the thermal part of the stress describing the 
physics of the glass transition, while aj is the athcrmal 
part of the stress accounting for the jamming transition. 
Finally, rj s is the solvent viscosity, so that r] s j is the stress 
stemming from the background solvent. This only be- 
comes relevant in the dilute limit but serves as a useful 
reference value when considering experimental data. 

Our goal with the model in Eq. © is not to make 
new types of predictions for the functional form of the 
flow curves for systems near glass and jamming transi- 
tions, but rather to explore the interplay between both 
physics. Thus, we make use of previous theoretical work 
and choose simple functional forms that most efficiently 
capture the physical features of the flow curves associated 
with each of these transitions. 



2. The glass contribution 

First, we specify the thermal part of the stress, ctq. 
We write it as the sum of the zero and the finite shear 
rate parts: 



(9) 



(jTTG(tp))- 1 + (1 + pgUtt)^)- 1 

Although compact, this expression incorporates quite a 
number of physical features of the rheology of glass- 
forming suspensions. In this equation, Yq and pq are 
two dimensionless numerical prefactors, with no deep 
information content. This expression is constructed to 
produce flow curves with no yield stress (defined as the 
7—^0 limit value of the shear stress) at low density, and 
a finite value for the glass phase when density is larger 
than the glass transition density, tpQ, which is assumed 
to be sharply defined. This last point is discussed further 
in Sections III CI and III El below. 

We note to start with that the shear rate and shear 
stress in Eq. §§§ are expressed in units appropriate for 



Brownian suspensions, namely <tt and Tr> It is mainly 
via these elementary units tt and <tt that the temper- 
ature dependence of the flow curves enters (apart from 
a minor temperature dependence of the transition den- 
sity, see Sec. Ill Cl below). The microscopic expressions of 
these units are given by Eqs. ([31 [5]) for a system obeying 
Langevin dynamics, as in our numerical simulations. In 
experiments, one might want to replace the damping co- 
efficient £ in this expression using instead Stokes' law to 
obtain the microscopic time scale: 



TT 



3irr] s a 3 



(10) 



Having set the correct scales, we now describe the 
functional form of the flow curves predicted by Eq. ^ . 
For densities below the glass transition, tp < ipc, the 
yield stress is zero, and we impose accordingly oqy{p < 
Pg) = 0. Thus, only the second term on the right 
hand side needs to be discussed. For small shear rates, 
7 <C (ttG(v?)) , the first term in the denominator dom- 
inates and the equation becomes ctq ~ YG<JTTrG((p) A f. 
This represents Newtonian behavior with the thermal vis- 
cosity rfr(ip) given by 



rft{<p) = 3iiYGG(tp)r) s 



(11) 



This expression shows that the dimensionless function 
G(tp) controls the rapid growth of the viscosity on ap- 
proaching the glass transition, and thus also corresponds 
to the growing equilibrium relaxation time of the un- 
sheared Brownian suspension at density tp, via the rela- 
tion T a RJ G((/?)tt. 

At larger shear rates, the flow curve © becomes 
o~g ~ Ygctt(1 + Pg(7 t "t) qg ); here we assume that the 
exponent qq takes a value below unity. This expression 
describes the strong shear-thinning regime obtained in 
dense fluids when the shear rate competes with the slow 
glassy dynamics and drives the system away from equilib- 
rium. This competition is modelled by a plateau regime, 
ctg ~ Yq <jt, when 7 -C 1 , finally followed by a different 
shear-thinning behaviour, <tg ~ ^ / gPG( j t(7''t) qg . These 
two shear-thinning regimes result from the well-known 
existence of two separated relaxation processes (a and ft 
relaxations) in highly viscous fluids [27], 

In the glass phase, tp > tpa, we impose an infinite 
value of the Newtonian viscosity using G(tp) = 00, such 
that the first term in the denominator disappears from 
Eq. © and the flow curve simplifies to ctg/ctt = [Yq + 
&gy(<p)] +YGPG(i T T) aa ■ This describes the existence of 
a finite yield stress, 



(12) 



whose scale is set by the thermal stress. At finite shear 
rates the rheology is of Herschel-Bulkley type [l|, with 
shear-thinning exponent «g- 

The resulting functional forms for the glass contribu- 
tion are shown graphically in the bottom left sector of 
Fig. [21 which corresponds to the thermal sector of the 



5 



io lv 
io 8 



«5 io 4 

I 



10 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 




















//^ Jamming 
transition 












/ Glass 1 
/ transition! 

< i , i . i 





icr 10" 10 Z 10' 



FIG. 2: Flow curves predicted by the additive model at low 
temperatures, feT/e = 1CF 10 , when glass (bottom left) and 
jamming (top right) sectors are well separated. Notice the 
different functional forms of the flow curves near the criti- 
cal densities tpa and ipj (highlighted with dashed and dotted 
lines), reflecting the distinct microscopic dynamics associated 
with thermal and athermal situations. 



rheology predicted by the present model. We defer a dis- 
cussion of the volume fraction dependences of ctgy and 
G to Sec. EC] below. 



3. The jamming contribution 

We now turn to the jamming contribution to the shear 
stress in Eq. ((8|) , crj . We express this athermal part of 
the stress as: 



CM 0,7) 



CO 



( 7 r J(^))- 1 + (pj(jT ) a >y 



(13) 



In this expression, we have introduced two dimensionlcss 
parameters, Yj and pj. In contrast to the glass con- 
tribution, the stress and time units are now cr and To, 
whose microscopic expressions were given in Eqs. JI1[G]). 
Therefore, temperature does not enter the stress contri- 
bution (|13|) . as expected from the athermal nature of 
jamming physics. 

Expression (fTB")) has many similarities with Eq. © but 
is qualitatively slightly simpler, as we now describe. For 
densities below the jamming density, we impose ajy(ip < 
tpj) = 0, such that a Newtonian viscosity emerges in 
the low shear rate limit, 7 <C (r J (if))" 1 . Here crj w 
Yj<jqTq J (<fi), which defines the athermal viscosity 



(14) 



This equation shows that the dimensionlcss function J(ip) 
now controls the divergence of the Newtonian viscosity on 
appraoching ipj. The difference with the glass regime is 



seen at larger shear rates because the constant contribu- 
tion in the denominator of Eq. ([§]) is absent in Eq. f| 13[) . 
As a result, Eq. (|13|) shows a simple power law shear- 
thinning behavior crj w Yjo~opj( A /To) aj , instead of the 
plateau in Eq. (|9]). This reflects the absence of distinct 
a and (3 relaxations in athermal systems [22j . 

In the jammed phase, ip > ipj, we impose a nonzero 
yield stress, ctjy > 0, and an infinite viscosity via J(tp) = 
00. As a result, we obtain (T j /do = cMy(^)+XiPj(7'?o) Q!j ) 
which describes the existence of a finite yield stress, 



yield 



(15) 



Its units are given by the athermal stress. Again at finite 
shear rates we have Herschel-Bulkley rheology, with the 
shear-thinning exponent now aj. 

The resulting functional forms for the jamming contri- 
bution are shown graphically in the top right quadrant of 
Fig. [21 which corresponds to the athermal sector of the 
rheology predicted by the present model. 



C. Density dependence: Choice of fitting functions 

In both thermal and athermal expressions for the shear 
stress, we have described a change from fluid to solid 
behaviour signalled by rapidly growing viscosities con- 
trolled by the dimensionless functions G(ip) and J(tp), 
and emerging yield stresses controlled by the functions 
cgy(v) an d ojY(y). Within our model, these functions 
entirely control the density dependence of the flow curves. 
In this subsection, we describe the physics embodied by 
these quantities, and motivate the specific choices we 
made to fit numerical and experimental flow curves in 
the present paper. 



1. Diverging Newtonian viscosities 

For thermal systems, the function G((p) controlling the 
growth of the Newtonian viscosity has been carefully an- 
alyzed in a large number of studies, as it quantifies the 
dynamic slowing down of a hard sphere fluid on its ap- 
proach to the colloidal glass transition 0, |28r]3fj| . As 
is well known from the literature of glass-forming ma- 
terials 0], this timescale growth cannot be represented 
by a single functional form because the slowing down 
of the dynamics crosses over from power law timescale 
increase at moderate volume fraction, to a steeper 'ac- 
tivated' exponential growth at larger density |30fl . This 
crossover can be understood theoretically as a change of 
relaxation mechanism from collective but nonactivatcd 
dynamics described by mode-coupling theory [3lj , to cor- 
related activated events closer to the glass transition 0] . 

In practice, however, given the limited range of viscos- 
ity data that is accessible in typical numerical and ex- 
perimental work, a single functional form is often a good 
enough approximation to describe the data. Therefore, 



G 



in the following we shall describe the data using a power 
law divergence: 

G(<p) = h G (vG - v)~ 7G , (16) 

where Jiq is a numerical prcfactor. The location of the 
glass transition, tpa, and the associated exponent, 7g, 
are fit parameters that serve to describe the increase with 
density of the thermal viscosity rjx((p). As is clear from 
the above description, the value extracted for the loca- 
tion of the glass transition has the same physical con- 
tent as the location of the mode-coupling temperature 
for standard glass-forming materials; it should not be 
confused with the location of a genuine dynamic or ther- 
modynamic singularity [32j . 

On the other hand, much less is known about the den- 
sity dependence of the athermal viscosity t]q, controlled 
in our model by the dimensionless function J(f). Theo- 
retical [H, HU and experimental [H, HH studies seem to 
indicate that a power law divergence describes the results 
well, and accordingly we use 

J{<p) = h 3 {<p 3 -<p)-n, (17) 

with hj a numerical prefactor. The exponent 7j is found 
numerically and experimentally to be close to 7,1 = 2. 
The critical density, if] , appearing in the expression (|17[) 
represents the location of the jamming transition. 

Although Eqs. ((TBI [T7|) appear similar, we note that 
the former has been derived from a microscopic perspec- 
tive HJ , while the second one is only supported by exper- 
imental data or detailed theoretical and numerical anal- 
ysis performed over a limited range of densities 0, [34j . 
Therefore, while empirically robust, the status of the al- 
gebraic divergence for non-Brownian hard sphere suspen- 
sions remains to be clarified. A second difference between 
the two expressions concerns the nature of the critical 
densities. While we described ipc as the crossover den- 
sity emerging from a mode-coupling theory analysis of 
the data, an actual divergence is expected in the ather- 
mal limit, such that <pj is not simply a crossover but 
a genuine critical density Q. However, its location re- 
mains ambiguous: the nonequilibrium nature of the jam- 
ming transition implies that it depends on the specific 
details of the procedure used to bring a system close to 
jamming [35) . This implies, for instance, that the values 
for (pj obtained by slow or fast compressions of random 
sphere packings do not coincide with the critical den- 
sity measured under shear [36j . Thus, ipj in Eq. (I17|) is 
very specifically defined as the critical density that con- 
trols the divergence of the Newtonian viscosity of non- 
Brownian particles in the hard sphere limit |l6l . Il9j . 

2. Emerging yield stress 

We now discuss the behaviour of the yield stress for 
glass and jammed phases, which is also characterized by 
a number of singular behaviours, as demonstrated by the 
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FIG. 3: Yield stress of harmonic spheres, measured numer- 
ically at fcsT/e = 10~ 4 , 1CT 5 , 1(T 6 and 1CT 7 (from top to 
bottom) . These results are extracted from the 7 — >■ limit of 
the flow curves in Fig.[TJ The red lines are fits using Eqs. (|21l 
1221) . The three lines for T = 10~ 7 represent fits to the sim- 
pler functional forms in Eq. (fl8|) [(i), dashed], Eq. (TT9)) [(ii), 
dash-dotted], and Eq. ([2(2) [(iii), dotted]. 

numerical measurements performed on harmonic spheres 
at various temperatures and reported in Fig. |3] These 
yield stress data are obtained from the zero shear rate 
extrapolation of the flow curves shown in Fig. [TJ 

As can be seen from Fig. 02 the yield stress data show 
three distinct characteristic behaviours, which are more 
or less pronounced depending on temperature, and thus 
one can use various levels of sophistication in the descrip- 
tion of these data. 

The first observation is related to the appearance and 
growth with density of a yield stress above the glass tran- 
sition lpq . Within our model, this behaviour is controlled 
by the term <jgy(<p), see Eq. (|12[). A possible functional 
form is 

<Jgy{¥ > </>g) = A g {lp - tp G f G , (18) 

where Aq is a prefactor and /3q an exponent which char- 
acterizes the growth of the yield stress with density. This 
form is consistent with the description provided by the 
extension of the mode-coupling theory to sheared sus- 
pensions [l4[, where the exponent is given by (3q = h- 
The data in Fig. [3] indicate that Eq. (fTS)) is sufficient to 
describe the yield stress of harmonic spheres when the 
temperature is high enough, e.g. fc^T = 1CT 4 , so that 
the jamming transition does not play any role. Note, 
however, that a precise determination of the exponent in 
Eq. (|18|) is not possible, as it is difficult to obtain un- 
ambiguous yield stress data very near ipc due to the 
crossover nature of the mode-coupling transition: the 
flow curve at the fitted Lpa would eventually become New- 
tonian at low enough shear rates. 

By contrast, the jamming density becomes relevant 
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when temperatures becomes much smaller, and this af- 
fects the yield stress on both sides of the jamming transi- 
tion. For tp < ip -f and T — > the behaviour of hard sphere 
glasses is recovered, for which the yield stress diverges as 
tpj is approached from below 37J. One can model this 
behaviour using again an algebraic dependence, 



CTGY^ < V>j) = B G (lfij - tp) 



-Pgj 



(19) 



where Bq is a numerical prcfactor. We are not aware of 
specific theoretical predictions for the value of the expo- 
nent /?gj governing the yield stress divergence, but one 
might expect that it is the same exponent which also 
controls the divergence of the pressure in compressed 
hard spheres, for which the value = 1 is well doc- 
umented HI HI, HI]. 

The third characteristic density dependence is ob- 
tained above the jamming transition in the T — > limit. 
In our model this is entirely controlled by the term a jy; 
see Eq. (fT5|) . The emergence of a yield stress in athermal 
soft sphere packings has been described in previous nu- 
merical work as a continuous power law growth [l|| H3] , 



(20) 



Note that there is no need to introduce a new prefactor 
in this expression, as the second term of Eq. (|13|) already 
contains two adjustable prcfactors in the denominator, 
pj and hj from J (tp). The expression (|20[) is known to 
be sufficient when describing fully athermal assemblies of 
soft particles [l6|. The exponent /3j has been found to 
be close to (3j = 1, although small deviations from this 
simple value of /3j have also been discussed . 

We have represented the three fitting functions in 
Eqs. (fT8l HH [20]) as lines going through the numerical 
results in Fig. [3J demonstrating that each of these ex- 
pressions can describe a limited range of densities fairly 
well; each therefore correctly captures a different physical 
regime of the yield stress behaviour of harmonic spheres. 

However, as is clear from Fig. [3j for finite but very 
low temperatures all three regimes affect the density de- 
pendence of the yield stress. Therefore, to describe the 
functional form of such low-T data correctly one must 
introduce more complicated model equations that inter- 
polate smoothly between the various regimes. A possible 
expression is as follows: 



0"GY 



ojy 




l (p-p G )-^(pj~ip)^ 



+ Y J (tp-tp J )^. 



(22) 



In these equations we have used the product Y G (tp — 
PgY g {pj—p)~P G} in the thermal part to interpolate con- 
tinuously between expressions (|18|) and (|T9| in between 
the glass and jamming transition densities. This has the 
added benefit of reducing the number of adjustable nu- 
merical prcfactors by one, by effectively replacing Aq 



and Bq with Y G . On the other hand, we have also in- 
troduced the constant Yqj in order to extend the T = 
singularities in Eqs. (fT9| and ([20]) to nonzero tempera- 
tures while connecting them smoothly. Indeed, one veri- 
fies easily that exactly at tpj the expressions (|21I22[) give 
the same result a- ylcld = ogyot = cjYfo = Igj V TT(Jo ' 
Therefore, the dimcnsionless parameter fc^T/e is again 
the key parameter controlling the typical stress scale at 
the crossover between thermal and athermal regimes, via 
y^Jr/fTo = y/ksT/e. Note that the complex pattern of 
these yield stress data near jamming resembles the be- 
haviour found for the pressure in thermalized packings 
of soft spheres near the jamming transition [40l |. 

For definiteness we note that in writing Eqs. (|2"T1 1221) . 
we have assumed that ojy is zero for tp < tpj while cqy 
vanishes for tp > tp j . The sudden drop of ctgy to zero at 
tpj does not, of course, have any physical meaning; only 
the sum (jt&gy + cocjy matters in our model and this 
is a smooth function of the volume fraction tp. 

In Fig. [3J the full lines through the numerical data are 
obtained by using simultaneously Eqs. ([2T1 [22| . which 
clearly gives very satisfying results. The fitting parame- 
ters used to describe the simulations are summarized in 
Table [I] below. We emphasize that the apparent complex- 
ity of these expressions and the relatively large number 
of adjustable parameters are needed because we want to 
capture in a single set of equations an unusually large 
number of physical phenomena pertaining to the physics 
of both glass and jamming transitions, for both thermal 
and athermal systems, soft and hard particles. 



D. The additive model at work 

Having carefully analyzed the behaviour of the yield 
stress in the previous section, we are now in a position 
to combine all the elements of the additive model to de- 
scribe the flow curves obtained numerically for harmonic 
spheres at various temperatures and densities. 

To this end, we add the glass and jamming contribu- 
tions, Eqs. (0 [T3j) respectively, to the shear stress @; 
the final solvent contribution can be safely neglected at 
the volume fractions of interest. Into these expressions 
we insert the fitting forms for the viscosity divergences, 
Eqs. ([TBI [TT]) , and the emerging yield stresses, Eqs. ([2"T1 
23) • The details of the chosen parameter values are given 
in Table [TJ below. 

This procedure yields the flow curves shown as continu- 
ous lines in Fig.[JJ In practice, we first tuned the parame- 
ters to reproduce the simulation results at fc^T/e = 10~ 6 . 
We then used these parameters to fit all other temper- 
atures. The figure demonstrates that the complex be- 
haviour observed in both glass and jamming limits, as 
well as the crossover between the two regimes, are well 
described by the additive model. This supports the va- 
lidity of our analysis in terms of a linear combination 
of independent contributions respectively stemming from 
glass and jamming physics. 
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E. Beyond the simplified mode-coupling 
description 

We have commented already several times on the fact 
that our additive model uses a simplified representation 
of the glass transition and its effect on shear rheology, 
which is in the spirit of mode coupling theory. Before 
proceeding to the analysis of experimental data, we want 
to clarify that this simplification is not conceptually re- 
quired by the additive model, and could be removed if 
desired. 

At first sight our way of writing the glass contribution 
to the stress, Eq. ©, seems to be tied to the existence of 
a specific volume fraction tpQ above which the yield stress 
cgy(v?) becomes nonzero. However, this contribution can 
be written in an equivalent form as 

a G (T,ip,j) = cr GY {f) +Y G 

(Tt ' (7TTGM)- 1 + (1 +PG(^)(7TT) aG ) ^ 

Here we have allowed a density-dependence of the nu- 
merical coefficient p G (ip). If this is chosen as Pg^) = 
Pg^g/(1g -I- ctgy (<£>))> it is a simple matter to check that 
© and ([25)) arc identical. This is because for <p < ipc we 
have ctgy(v) = 0, while for tp > tpQ, G(tp) is infinite and 
so the first term in the denominator vanishes. 

The form of the glass stress is now easy to gen- 
eralize to more sophisticated representations of the be- 
haviour of the glass transition. One could represent the 
crossover to activated relaxation processes by keeping 
G((p) finite but making it cross over to an exponentially 
fast increase above ipc- In the same vein, ogy(¥>) could 
be given a smooth, non-singular onset around ipc, and 
one could also allow a more general density dependence 
for pg (</?)• Conceptually, these generalizations are im- 
portant: if G(<p) remains finite, the model never predicts 
a genuine yield stress for tp < ipj and there is always a 
finite Newtonian viscosity, given by ([TT]) . A final possi- 
bility would be to incorporate a genuine glass divergence 
located at the 'ideal' glass transition density \§L |32|. 

For representing real flow curves, on the other hand, 
it is essentially irrelevant whether G(cp) has a true di- 
vergence or not. In the latter case, the value G(ip) will 
still become so large near <pc that 7 3> (ttG^))^ 1 for 
all accessible shear rates. The predicted flow curve then 
exhibits an effective yield stress, given by (fT2|) as before. 
Our simpler representation of G(tp) avoids introducing 
additional parameters for the divergence above ipc , which 
cannot be determined with any accuracy from the avail- 
able rheological data and would not improve the quality 
of our fits. 



III. ANALYSIS OF EXPERIMENTAL RESULTS 

A. Overview 

In this section, we will apply the set of equations de- 
scribed in Sec.[H]to analyze experimental data stemming 



from a number of different sources. Our main goal is to 
demonstrate how to disentangle glass and jamming rhe- 
ologics in experimental work, and detect whether there 
exist experimental systems where the crossover between 
both sectors revealed by our numerical work is relevant. 
The study of such systems would be useful in assessing 
the validity of the key aspect of the model — the additiv- 
ity hypothesis. Another expected outcome is a clearer 
understanding and classification of the nature of the for- 
mation of amorphous solids in soft materials. 

For all systems considered, we shall need to determine 
the particle softness, the particle size and the solvent vis- 
cosity, which will serve to construct the elementary time 
and stress scales. In our additive model, the particle soft- 
ness ksT/e is the key parameter governing the structure 
of the flow curves. As mentioned above, very soft parti- 
cles will not be very sensitive to the jamming singularity. 
This is because the latter requires that the notion of par- 
ticle contact is physically meaningful, which is not the 
case when large thermal fluctuations are present. At low 
temperatures, on the other hand, the two transitions are 
well-resolved and belong to distinct stress and time sec- 
tors as demonstrated by the flow curves in Fig. [5J For 
fixed low particle softness, the role of the particle size and 
solvent viscosity is thus to determine which sector of the 
flow curves is actually observed in any given experiment. 

Before going into the details of the various systems, it 
is instructive to consider some numbers. In a typical soft 
matter experimental setup [l|, measurable shear stress 
and shear rate scales are a w 1 Pa and 7 w Is" 1 . For a 
particle diameter of the order of a = 0.1 /*m and solvent 
viscosity r] s = 1 mPa • s (this is the typical value for wa- 
ter), one obtains ct/ctt = era 3 /ksT as 1 and 7TT ~ 10~ 3 . 
This implies that the thermal sector will be explored ex- 
perimentally, and glass transition physics should be dom- 
inant. 

On the other hand, when the particle diameter is 
a = 10 /im with the same solvent, otherwise identical ex- 
perimental conditions will correspond to era 3 /ksT ss 10 6 
and 7TT ~ 10 3 . Such an experiment would probe the 
athermal sector of the flow curves, where a jamming tran- 
sition can be observed. 

This little exercise shows that an interesting inter- 
play between glass and jamming physics can be observed 
only in an experiment with particle sizes of the order of 
a ss 1 ^tm, and for particles which are not too soft. Wc 
shall see below that emulsions seem to represent the best 
compromise. In the next subsections, we follow the above 
approach: we shall redraw several available experimental 
flow curves using, for definitcness, units appropriate for 
thermal systems; we then fit these flow curves using the 
model equations described in Sec. QTJ 



B. Glass physics: PMMA colloids 

In this subsection, we analyze the flow curves reported 
in Ref. j4l[ for nearly hard sphere colloids. This PMMA 
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FIG. 4: Experimental flow curves of a nearly hard sphere 
PMMA colloidal suspension, taken from Petekidis et al. [4l| . 
are shown using thermal units. The full lines are the fits us- 
ing our additive model using fcsT/e = 10 -8 and ipG = 0.575; 
see also Table [I] For the first glass flow curve, we show the 
decomposition of the stress into its glass (dotted) and jam- 
ming (dashed) contributions. The thin dashed line represents 
the dilute limit. Clearly, the fluid to solid transition is a col- 
loidal glass transition, with little influence from the jamming 
physics. 



colloidal suspension is a very important experimental sys- 
tem whose rhcology has been studied in much detail 
This is because PMMA colloids are considered as a ver 
good experimental realization of the hard sphere fluid 
which is itself an important model system for the statisti- 
cal mechanics of simple fluids [42j. In particular, PMMA 
colloids have been used extensively in experimental stud- 
ies of the hard sphere glass transition [f| I28U301 |43| . In 



the specific study in Ref. [41[, the particle diameter is 
a = 0.36 /im with a size polydispersity of about 12 % to 
prevent crystallization. We can thus estimate the ther- 
mal stress scale as ctt = fcsT/a 3 = 0.0825 Pa. Since the 
free diffusion constant D of this colloid is also reported 
in Ref. [4l[, the thermal time scale for this system can 
be estimated by tt = a 2 /Do without using the Stokes' 
law Eq. (flU)) . The estimated value is tt = 0.158 s. We 
show the measured flow curves using these thermal units 
in Fig.|4j This representation quickly establishes that the 
flow curves are essentially located in the thermal sector 
described by our glass/jamming rheology model. In par- 
ticular, the transition from a Newtonian fluid to a yield 
stress solid occurs at low Peclet number. 

We have fitted the flow curves to the additive model. 
The fitting parameters are summarized in Table HI Note 
that the particle softness /csT/e = 10~ 8 used in our anal- 
ysis is somewhat arbitrary, since this parameter does not 
strongly affect the flow curves in this thermal regime, 
as expected from the nearly hard sphere nature of the 
particles. As shown in Fig. 21 the transition between 



fluid and solid states is very well described by the glass 
rheology, and the fits indicate that the transition occurs 
for the volume fraction ipQ = 0.575. This is consistent 
with the analysis performed in the original paper [4lj . 
This value is also consistent with a mode-coupling analy- 
sis of the microscopic dynamics of colloidal hard spheres, 
for instance using light-scattering [29l [30| or microscopy 
techniques [43| . Thus our analysis confirms that the dy- 
namic range probed in a typical rheology experiment is 
not broad enough for the crossover to activated dynam- 
ics at large density (which has been observed numeri- 
cally [HI and using other experimental approaches (30[) 
to become apparent. 

Although it is clearly glass physics that controls the 
overall features of the flow curves, it is interesting to ask 
whether jamming physics also manifests itself. We argue 
that in Fig. @] there are two distinct aspects where this 
is the case. First, we note that the yield stress reaches 
the value a yield /<7T ~ 10 2 for the largest volume fraction, 
ip = 0.62. This large value implies that the density is not 
very far from the jamming density so that the divergence 
in Eq. (|19[) starts to become relevant and a simplified 
expression for the yield stress, such as Eq. (|T5|) . would 
not account for the experimental data taken deep in the 
glass phase. 

Accordingly, if the density is close enough to tpj , then 
the athermal Newtonian viscosity rjo(ip) should also start 
to become large, as it is controlled by a similar diverging 
expression, see Eq. ([TT|). From the above discussion of 
the model, see for instance Fig. [2 this viscosity is con- 
tributing to the flow curves at large Peclet number. This 
is demonstrated in Fig. U by the dashed line represent- 
ing the jamming contribution to the total shear stress, 
cj ~ ?7o(</?)7- However, while our model predicts that 770 
grows rapidly when ip — > ipj, such a growth is not ob- 
served experimentally, and there are clear deviations be- 
tween the fit and the data for the largest density in Fig.|U 
This might indicate that the PMMA colloids cease to be- 
have as nearly hard spheres at large densities and large 
Peclet numbers. 

We suggest that it would be interesting to repeat 
such steady state rhcological measurements with slightly 
larger PMMA particles so that large Peclet numbers and 
large densities are more easily studied. The crossover be- 
tween glass and jamming rheology should then become 
more apparent and its features could be elucidated ex- 
perimentally in this well-studied colloidal system. 



C. Jamming physics: Aqueous foam 

We now analyze the flow curves of aqueous foam re- 
ported by Herzhaft et al. [Hj]. Foams are considered as 
prototypical materials displaying a jamming transition, 
because their are typically made of non-colloidal soft bub- 
bles. It is worth mcntionning that the harmonic sphere 
model considered numerically in Sec. Ill Al was first de- 
vised to study the jamming rhcology of wet foams 24 1. 
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FIG. 5: The experimental flow curves of aqueous foam in 
thermal units, after Herzhaft et al. [ii] . The red curves are 
fit by the additive model with ksT/e = 2.10 -11 and ipj = 
0.55; see also Table U The fluid to solid transition occurs in 
the jamming sector, with negligible influence from thermal 
fluctuations. 



FIG. 6: Experimental flow curves of oil-in-water emulsion in 
thermal units, after Mason et al. [47j. The red curves are 
fits to the additive model with ksT/e = 3 ■ 10 -7 , ipG = 0.579 
and <pj — 0.62; see also Table[U While a glass transition is ob- 
served at low shear rates, the sharp increase of the yield stress 
in the glass phase is a signature of the jamming transition. 



The experimental system is composed of nitrogen bub- 
bles dispersed in an aqueous polymer solution. Flow 
curve were measured at various densities. However, the 
particle size also changes with the density. We extracted 
the average particle diameter for each density from the 
reported particle size distributions. Note that the typical 
diameter is rather large at 27 /mi, with some size polydis- 
persity. The reported solvent viscosity is rj a = 13 mPa-s. 
Using these values, we have estimated the thermal time 
scale and stress scale for each density. Typical values 
are tt = 5 • 10 5 s and <tt = 2 • 10~ 7 Pa, which are clearly 
outside the experimentally accessible windows for typical 
time and stress measurements. We represent the exper- 
imental flow curves using these thermal units in Fig. [5] 
As expected this produces large dimensionless numbers 
(such as 7Tt ~ 10 7 — 10 9 ). Clearly, these flow curves 
belong to the athermal sector, and should mainly be con- 
trolled by the jamming contribution to the shear stress. 

In the case of a bubble or droplet, we can estimate 
the particle softness using the surface tension A. The 
pressure inside a bubble is larger by AP = 4 A/ a than 
on the outside, AP being the Laplace pressure. This 
pressure difference acts as a repulsive force between two 
overlapping particles. When the overlap length is d, the 
interface area between the 'overlapping' bubbles is given 
by ira 2 d/2, and the repulsive force becomes 2irAad. In- 
terpreting this force as the derivative of a pair interaction 
potential, we estimate the softness as e = 2-kAcl 2 . This 
expression and the reported surface tension A lead us to 
estimate the dimensionless particle softness for this sys- 
tem to be ksT/e = 3 • 10 _u , which is indeed very close 
to the athermal limit in which jamming physics should 
be observed. 



We fitted the experimental flow curves to our model 
equation using this particle softness. The result and fit- 
ting parameters are shown in Fig. [5] and Table HI The 
transition is well fitted as a jamming transition with 
ipj = 0.55. Note that this value is much lower than the 
random close packing density usually quoted for spherical 
particles, ipj « 0.64. The reason of this deviation is not 
very clear. We can invoke the fact that the interaction 
between real bubbles is more complicated than in simple 
models of repulsive spheres, or the idea that additional 
microscopic dissipation channels exist in the real mate- 
rial, due for instance to some internal degrees of freedom 
of the bubbles. We note, however, that in other studies 
the jamming transition of a real foam was identified and 
located very close to the jamming transition occurring 
in simple harmonic sphere models, both at the level of 
static [45| and rheological [46| properties. 



D. Exploring the crossover: Oil-in-water emulsions 

We next focus on the flow curves obtained for oil-in- 
water emulsions, reported by Mason et at [47j . The sys- 
tem is composed of oil droplets stabilized with sodium 
dodccyl sulfate, which are dispersed in water. Although 
experiments with different particle diameters were per- 
formed [471 ] , we only analyze the results for droplet diam- 
eter a = 0.5 /mi (with less than 10 % size dispersity), be- 
cause the most complete set of flow curves were reported 
for this specific particle size [47| ■ Using the water viscos- 
ity value of t) s = 1 mPa-s, we estimate the thermal time 
and stress scales to be tt = 0.25 s and or = 3 • 10~ 2 Pa, 
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respectively. These values are well within the measur- 
able range and then we expect aspects of the glass tran- 
sition to be relevant, although large Peclet numbers can 
certainly be accessed too. However, it is important to 
note that, in contrast to PMMA colloids as discussed in 
Sec. IIII Bl emulsions are made of soft droplets which can 
therefore easily be compressed above the jamming (or 
'random close packing') density. This implies that as- 
pects of the jamming transition are also potentially im- 
portant for this system. The fact that features of both 
transitions are relevant for these emulsions is already ap- 
parent from the original experimental papers: the lin- 
ear rheology of these systems was interpreted using the 
concepts of a and /3 relaxations and power law scalings 
near the glass transition density ipQ inspired by mode- 
coupling theory 48[ , while the density dependence of the 
yield stress was fitted to a power law near the jamming 
density tpj in a separate article (47} ■ Our additive model 
represents an ideal framework to unify and rationalize 
these findings. 

The complex rheological features summarized above 
can be seen in Fig.[6j where the experimental flow curves 
are shown in thermal units. The system behaves as a 
Newtonian fluid at low enough density, tp < 0.57, and 
as a yield stress fluid for larger densities. However, the 
behaviour of the yield stress in the glass phase is clearly 
nontrivial, as it increases sharply with tp in the range 
tp ss 0.60 — 0.62 to reach large dimensionless values, 
a y' ,cld /a T « 4 • 10 2 for tp = 0.65. This final shear stress 
value should belong to the athermal regime, possibly sug- 
gesting a crossover between the glass and jammed phases. 

To confirm these qualitative conclusions, we fit the 
experimental flow curves using the additive rheological 
model. To this end, we first estimate the particle softness 
from the surface tension reported in the experimental ar- 
ticle (47| , using the method outlined in the previous sub- 
section for aqueous foams. We obtain ksT/e = 3 • 10~ 7 . 
Interestingly, this value is much larger than the one ob- 
tained for both foams and PMMA colloids, but remains 
at the lower end of the range of temperatures simulated 
in our numerical simulations [22J; see Fig. [TJ We then 
fitted the experimental flow curves to the model equa- 
tion using this softness value, and obtained the fits shown 
with lines in Fig. El The corresponding fitting parameters 
are listed in Table. |U We obtain the glass and jamming 
transition densities at tpQ = 0.579 and pj = 0.62, these 
values b eing consistent with previous experimental anal- 
ysis [13, l48l ] We emphasize that both the fluid to solid 
transition observed at low Peclet numbers, and the sharp 
increase with density of the yield stress arc well fitted by 
our model. This allows us to conclude unambiguously 
that the transition which appears near ip = 0.58 is a col- 
loidal glass transition, while that the sharp stress increase 
at higher density near tp ~ 0.62 is a clear signature of the 
jamming transition. We stress that to obtain a quan- 
titative description of these flow curves, it is absolutely 
necessary that both glass and jamming contributions to 
the shear stress are combined in a unified model as in 
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FIG. 7: Viscosity vs. shear flow curves measured (a) exper- 
imentally for oil-in-water emulsions [47| and (b) numerically 
for harmonic spheres at ksT/e — 10 -6 . We use similar di- 
mensionless scales to emphasize the similarity between the 
two situations, suggesting that emulsions are an ideal experi- 
mental system to study the glass-jamming crossover. 



Eq. dSJ. Therefore, the excellent fits in Fig. |H] provide 
strong experimental support for the additive model and 
theoretical analysis offered in the present work. As a di- 
rect physical corollary, we also conclude that the shear 
rheology of micron-size emulsions is quite complex, as the 
flow curves reveal all the features described by the most 
complex version of our additive rheological model. 

As shown in Table |U the parameters used to describe 
the emulsion data are very similar to the ones used to 
fit the flow curves of harmonic spheres in Fig. [T] This 
implies that the two sets of data should in fact be quan- 
titatively very close. To show this more directly, we re- 
plot the flow curves for both these systems in a slightly 
different representation in Fig. where the shear rate 
dependence of the effective viscosity, 77(7) = cr( A f)/ A f, is 
shown. We make the experimental viscosity dimension- 
less using the solvent viscosity rj Sl and combine Stokes' 
law in Eq. (TTU)) with Eq. © to obtain the bare vis- 
cosity rj s = £/(37ra) in the numerics. The results in 
Fig. [7] arc clearly very similar. This means that, some- 
what surprisingly, modelling of real emulsions as har- 
monic spheres with overdamped Langevin dynamics is 
quantitatively very accurate. Our results also suggest 
that an interesting athermal Newtonian behavior char- 
acterized by the viscosity 770 would appear experimen- 
tally if higher shear rates could be studied. (Note that, 
unusually, the range of shear rates studied is actually 
broader here in the numerical work than in the exper- 
iments). This second viscosity should be directly con- 
trolled by the jamming transition, and emulsions could 
thus be used to improve our understanding of its density 
dependence, which is usually studied using non-Brownian 
granular suspensions [H,[3J)- We suggest that a more ex- 
tensive rheological investigation of emulsions, perhaps by 
varying the droplet size in a more systematic manner to 
gradually access larger Peclet numbers, would be very in- 
teresting. We have found literature data for much larger 
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particle diameters [49|, l5C§ , but the density range covered 
in these flow curves was too limited to perform a detailed 
analysis of their evolution using our model. 



E. PNIPAM microgels: Glass or jamming? 

In this final experimental subsection, we analyze the 
steady state flow curves obtained for two independent 
sets of similar PNIPAM microgels. These particles 
are currently the focus of a large number of investiga- 
tions [5lj ]. for at least two reasons. First, these parti- 
cles are made of microgels, and are therefore very soft. 
As such, they represent a new type of colloidal system, 
potentially very distinct from the more heavily studied 
hard sphere paradigm [52j. A second important feature 
is that PNIPAM microgels are highly sensitive to tem- 
perature, in the sense that a small temperature change 
induces a relatively large change of the particle size. This 
endows microgel suspensions with peculiar physical prop- 
erties [5lj . 

We have analyzed two sets of flow curves obtained 
experimentally. One is reported by Carrier et al. [53[ 
and was interpreted there in the framework of the mode- 
coupling theory extended to sheared fluids 1J], thus im- 
plying the assumption that the fluid to solid transition in 
microgel suspensions is a colloidal glass transition. The 
second one is reported by Nordstrom et al. [2l|. In this 
article, the flow curves are scaled using a procedure that 
was first employed in the context of the jamming tran- 
sition of soft spheres [l6( , and the authors interpret the 
fluid to solid transition as the result of jamming. Our 
model, based as it is on a careful distinction between 
glass and jamming physics, is thus a natural tool to re- 
solve the conflict between two opposing interpretations 
of the data. 

Like for other systems, our first task is to rcplot the 
flow curves in dimcnsionlcss units, using the particle size 
and solvent viscosity to determine the relevant thermal 
units; see Fig. [8] In the first experiment the au- 
thors had reported their flow curves using similar ther- 
mal units. However, the particle radius was used as a 
microscopic length scale. To be consistent with the anal- 
ysis performed elsewhere in the paper, we replottcd the 
results using the particle diameter as the unit length. In 
the second experiment [2l[, the results are reported in 
units of Pa and s _1 for shear stress and shear rate, re- 
spectively. The authors also reported the particle diam- 
eter as a function of temperature, and the corresponding 
packing fraction [63| . Using these diameters and the re- 
ported solvent viscosity, r) s = 1 mPa-s, we converted these 
flow curves to the dimensionless representation shown in 
Fig. [5J Once this rescaling is performed, it becomes ev- 
ident that the experimental results obtained by the two 
research groups extend over a very similar range of shear 
stresses and shear rates. Thus the simple dimensional 
procedure indicates that the physics probed in these two 
experiments is the same. In particular, we notice immc- 
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FIG. 8: Experimental flow curves of PNIPAM microgels of 
(left) Carrier et al. [k| and (right) Nordstrom et al. [2ll ]. 
using thermal units. The red curves are fit by the additive 
model with ksT/e = 10 -4 in both cases; see also Table|U The 
corresponding glass transition densities are <^>g = 0.68 (left) 
and ipG = 0.638 (right). 



diately that the experiments are performed over a range 
of shear rates corresponding to relatively small Peclet 
numbers, down to P c = 10~ 4 . This is because the par- 
ticle size is small enough for thermal fluctuations to be 
relevant over the typical timescale of the experiments. 
Note that when converted into dimensionless units using 
the athermal timescale To, as done in Ref. [2l|, the typ- 
ical time scales extracted from the flow curves become 
unphysically large, t/tq ~ 10 7 — 10 11 , again indicating 
that the experimental data do not belong to the athermal 
jamming sector. 

An independent confirmation of the relevance of ther- 
mal fluctuations for PNIPAM microgels is obtained from 
the linear viscoelastic experiments performed by Carrier 
et al. 53]. Their results show that the loss modulus 
G"(oj) becomes larger than the storage modulus G'(oj) 
at low frequency. This indicates that spontaneous relax- 
ation induced by thermal motion plays an important role 
in the rheology. Qualitatively similar results [48[ have 
been obtained, for instance, for the emulsion system dis- 
cussed in Sec. IIII Dl 

In order to fit the experimental data on PNIPAM mi- 
crogels to our additive model, we have to estimate the 
particle softness ksT/e. However, it is not obvious that 
the interaction between microgel particles can be approx- 
imated with a soft harmonic repulsion. Fortunately, a 
quantitative test of this assumption was performed in 
a detailed analysis of the vibrational properties of the 
amorphous solid phase of PNIPAM particles [5J| . For a 
similar system this study establishes the quality of the 
mapping and obtains ksT/e = 2 • 10~ 5 . A slightly larger 
value, near kgT/e = 10~ 4 , was later determined through 
a theoretical analysis of the mean-squared displacements 
of PNIPAM particles [H]. 

Using ksT/e = 1CP 4 for the value of the particle soft- 
ness, we then fitted both sets of flow curves to our ad- 
ditive model. The result of this procedure is shown in 
Fig. [51 The quality of the fits is similar to the one ob- 
tained within the mode-coupling theory framework [53| . 
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TABLE I: Summary of fitting parameters used in the paper. For emulsion, we used the same values of Yj and pj as in the 
original article [4?| . 



System 




pj 


ho 


hj 


OLG 


QJ 


7G 


7J 


Y G 


Y G 


Y G j 


Yj 


/3g 


/3gj 


At 


Simulation [22] 


7 


4 


0.03 


0.1 


0.3 


0.5 


2.2 


2.0 


0.15 


0.25 


0.02 


0.02 


0.6 


1.0 


1.0 


PMMA colloids [41] 


3 


4 


0.03 


0.3 


0.3 


0.5 


2.2 


2.0 


1.5 


6.0 


0.02 


0.02 


0.6 


1.0 


1.0 


Foam [44J 


7 


12 


0.03 


0.08 


0.3 


0.6 


2.2 


2.0 


0.15 


0.25 


0.02 


0.01 


0.6 


1.0 


1.0 


Emulsion [47] 


7 


0.5 


0.03 


0.06 


0.3 


0.5 


2.2 


2.0 


0.25 


0.65 


0.02 


0.016 


0.6 


1.0 


2.0 


PNIPAM [53] 


7 


1 


0.03 


0.07 


0.5 


0.4 


2.2 


2.0 


12 


30 


1.5 


0.03 


0.6 


1.0 


1.0 


PNIPAM [21] 


30 


1 


0.2 


0.005 


0.5 


0.4 


2.2 


2.0 


15 


500 


3.5 


3 


0.6 


1.0 


1.0 



As is clear also from the numerical flow curves obtained 
for a similar temperature in Fig.[TJ jamming physics plays 
virtually no role in these fits, and the main control pa- 
rameter to be adjusted is the glass transition density, 
which we estimate as ipp = 0.68 and tpQ = 0.638 for the 
first [53l | and second [2l| set of experiments, respectively. 
Therefore, we conclude that the transition between a vis- 
cous fluid and a soft solid observed at low shear rate in 
the data of Fig. [8] is a colloidal glass transition, thus fa- 
voring the interpretation of Carrier et al. [53j over the 
one of Nordstrom et al. [2l[ . 

We notice that another effect of the large particle soft- 
ness is that the glass transition density does not corre- 
spond to the one of the hard sphere system discussed in 
Sec. IIIIB1 but is considerably larger. This results from 
the fact that at finite temperatures microgcl particles 
can overlap while hard spheres cannot, thereby shifting 
the glass transition density to larger values, as observed 
numerically [23], HH and captured theoretically using a 
mode-coupling approach [56j. We emphasize that the 
quantitative similarity between the glass transition den- 
sity of harmonic spheres and the jamming transition of 
hard sphere particles is then nothing but a numerical co- 
incidence. 

Another difference between hard and soft particles lies 
in the value of the typical yield stress scales observed in 
the experiments, see Fig. [51 which is somewhat larger for 
the soft microgel particles. To account for this feature, 
we had to use a larger value of the parameter Yq (see 
Table U) which sets the scale of the yield stress in Eq. (f^Tj) . 
We note that a similar change of numerical prefactor by 
about an order of magnitude was introduced (with no 
discussion) in the mode-coupling analysis of Carrier et 
al. [11] . We speculate that the effect could be due to the 
influence of polymeric internal degrees of freedom of the 
microgcl particles. 

While the microgel rheology analyzed in this subsec- 
tion leads us to the conclusion that dense PNIPAM as- 
semblies should be considered as thermal glasses rather 
than athermal jammed solids, we suggest that performing 
experiments with larger PNIPAM particles would be ex- 
tremely rewarding, as this would produce an experimen- 
tal system where large Peclet numbers could be studied 
with particles that are very soft. This would allow the 
jamming transition and rheology of soft particle systems 
to be studied experimentally pa . l57l [58| . 



IV. SUMMARY AND PERSPECTIVES 

In this paper, we described in detail the 'additive' rheo- 
logical model first introduced in our numerical analysis of 
the interplay between glass and jamming rheologies [22| . 
and we used the model to revisit published experimental 
data. 

In our model, the stress contributions from thermal 
glass and athermal jamming physics are treated sepa- 
rately and we assume they can be linearly combined in 
a unified model. Because we want to describe the emer- 
gence of amorphous solids at large density for both col- 
loidal and non-Brownian repulsive particles, the main 
control parameters for the formation of soft solids are 
therefore the particle packing fraction, <p } and the parti- 
cle softness expressed in units of thermal energy, ksT /e. 
From the analysis of numerical data obtained for har- 
monic spheres at both finite and zero temperatures, we 
developed simple functional forms for the flow curves 
(J [pi) measured at a given state point, for the Newtonian 
viscosities t\t and 770 (in the fluid regimes), and for the 
yield stress cr ylold (in the solid phases). 

In particular, our additive model describes the emer- 
gence of solidity for soft repulsive systems as a function 
of temperature and volume fraction. As a result, we can 
use the yield stress fits shown in Fig. [3] to construct the 
three-dimensional 'jamming phase diagram' [ll[59| shown 
in Fig. [51 which represents the yield surface delimiting 
the solid and fluid phases for all values of the thermody- 
namic control parameters 22j. This surface shows that 
the emergence of solidity at finite temperatures is always 
controlled by the colloidal glass transition, the T = 
point being the only situation where the jamming tran- 
sition signals also the onset of solidity. 

The line describing the PNIPAM microgel data in 
Fig. M corresponds to high temperatures (soft particles), 
and therefore the jamming transition does not influence 
the rheology of these microgel suspensions. When tem- 
perature is decreased, a clearer signature of the jamming 
transition is observed as a sharp increase of the yield 
stress with density. In this regime, the jamming tran- 
sition appears as a change in the nature of the glass 
phase [io, [HJ , but it does not control the emergence of 
solidity which still occurs at the glass transition. This de- 
scription applies to emulsions, as shown in Fig.[SJ Finally, 
when thermal effects become negligible, the thermal yield 
stress might become too small to be detected experimen- 
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PNIPAM 




FIG. 9: Three-dimensional 'jamming phase diagram' showing 
the yield stress surface as a function of the thermodynamic 
parameters temperatures and density, in a dimensionless rep- 
resentation (particle softness ksT/e, volume fraction tp, and 
stress cra 3 /e). The three lines represent the location of the 
experimental systems discussed in Sec. Hm Foams are mainly 
sensitive to jamming physics, emulsions display an interest- 
ing interplay between glass and jamming transitions, while 
PNIPAM microgels undergo a colloidal glass transition. 




10~ 4 10~ 2 10° 10 2 10 4 10 6 10 1 



FIG. 10: Superposition of experimental flow curves for dif- 
ferent materials using thermal units. PMMA colloids with 
a = 0.36 fim, after Petekidis et al. [4l[. Foam with a — 27 fim, 
after Herzhaft et al. (44|. Emulsion(a) with a = 0.5 /jm, after 
Mason et al. PNIPAM(a) with a - 0.2 fim, after Carrier 
et al. j2l . PNIPAM(b) with a = 1.2 fim, after Nordstrom et 
al. [21 1 . Star polymers with a = 0.07 fim, after Koumakis et 
al. [62 1 . Emulsion(b) with a = 8pm, after Otsubo et al. [49l |. 



tally, and solidity genuinely emerges at the jamming tran- 
sition. This is the case for foams in Fig. [H] for which 
the glass 'wing' has negligible effects. Note that PMMA 
colloidal suspensions would appear at nearly the same 
temperature/softness as foams in the jamming phase di- 
agram of Fig. [§] However, with the particle size being 
much smaller than for foams, the yield stress emerging 
at the colloidal glass transition would easily be measured 
experimentally, and the measurements would stop as the 
jamming density is approached because the yield stress 
would seem to diverge there. 

As shown by the jamming phase diagram in Fig. [HI 
our analysis is useful in organizing the physics of differ- 
ent experimental systems. To confirm this, we have used 
our additive rheological model to analyze various exper- 
imental flow curves obtained for a variety of dense sus- 
pensions. The systems we focussed on were PMMA col 
loids [4l| , aqueous foam 



and PNIPAM microgels 



oil- in- water emulsions [47| . 
53|. We have also gathered 
experimental data from other sources, in particular ultra- 
soft particles composed of star polymers [62| , and data for 
emulsions with larger droplet sizes J4i|, but for brevity 
the results of our analysis have not been presented in 

secnni 

We showed that all the above experimental results can 
be successfully analyzed using the additive model. It is 
instructive to replot all data in a single figure using the 
dimensional procedure adopted throughout this paper, 
i.e. expressing stress and time scales in thermal units ctt 
and Tt, see Eqs. 0). These flow curves are collected 



in Fig. 1101 In this representation, the flow curves for 
PMMA colloids, star polymers, PNIPAM microgels lie in 
the same sector, which corresponds to the thermal sec- 
tor in our model; see Fig. [5] Therefore, the formation of 
amorphous solids in these systems stems from the physics 
of the colloidal glass transition. On the other hand, foams 
lie outside this regime and are controlled, accordingly, 
by the jamming transition. Interestingly, emulsions lie 
somewhat in between and so are influenced by both types 
of physics, as discussed in Sec. IIII Dl Note in particular 
that emulsions with larger droplet sizes, also shown in 
Fig. [TU1 could be useful systems to fill the gap between 
colloids and foams. While experimental studies of micro- 
gel particles have been interpreted from the point of view 
of the jamming transition [21], our analysis shows that 
for these soft colloidal particles the physics of jamming 
has, in fact, only a negligible effect. A similar conclu- 
sion has recently been reached based on the analysis of 
the short-time vibrational dynamics in the amorphous 
phase (55| . 

Our conclusion that glass and jamming rhcologics be- 
long to different sectors and contribute linearly to the 
shear stress is directly supported by the numerical flow 
curves obtained for harmonic spheres, and by the analysis 
of the oil-in-water emulsions in Sec. IIII Dl which clearly 
showed the complex features also observed in the simu- 
lations. We mentioned that similar indications are also 
found for PMMA colloids, in particular at large Peclet 
number and larger density, while microgel suspensions 
appear less well suited for a detailed experimental in- 
vestigations of the interplay between glass and jamming 
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rhcology. The data in Fig. [10] show that many current 
experimental systems are in fact too close to the thermal 
sector. Therefore, to fill the gap between small colloids 
and foams, measurements on sufficiently hard particles 
(PMMA colloids, emulsions) with particle sizes in the 
range a = 0.5 — 5/iin would appear ideally suited. We 
hope our work will stimulate experimental investigations 
in this direction. 
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